Random-roughness hydrodynamic boundary conditions 
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We report results of lattice Boltzmann simulations of a high-speed drainage of liquid films squeezed 
between a smooth sphere and a randomly rough plane. A significant decrease in the hydrodynamic 
resistance force as compared with that predicted for two smooth surfaces is observed. However, this 
force reduction does not represent slippage. The computed force is exactly the same as that between 
equivalent smooth surfaces obeying no-slip boundary conditions, but located at an intermediate 
position between peaks and valleys of asperities. The shift in hydrodynamic thickness is shown to 
depend on the height and density of roughness elements. Our results do not support some previous 
experimental conclusions on very large and shear-dependent boundary slip for similar systems. 

PACS numbers: 83.50. Rp,68.08.-p 



Introduction.— It has been recently well recognized 
that the famous no-slip boundary condition, for more 
than a hundred years applied to model experiments in 
fluid mechanics, reflected mostly a macroscopic charac- 
ter and insensitivity of old style experiments. Modern 
experiments concluded that although the no-slip postu- 
late is valid for molecularly smooth hydrophilic surfaces 
down to contact [iJ-EI, for many other systems it docs 
not apply when the size of a system is reduced to micro- 
and nano scales. The changes in hydrodynamic behavior 
are caused by an impact of interfacial phenomena, first 
of all hydrophobicity and roughness, on the flow. The 
effect of hydrophobicity on the flow past smooth surfaces 
is reasonably clear and suggests an amount of slippage 
described by the condition v s = bdv/dz where v s is the 
slip velocity at the wall, b the slip length, and the axis 
z is normal to the surface. The assumption is justified 
theoretically [4|-|7|and was confirmed by surface force 
apparatus (SFA) 4] , atomic force microscope (AFM) [l| , 
and fluorescence cross-correlation (FCS) Q experiments. 
Despite some remaining controversies in the data and 
amount of slip (cf. Q), a concept of hydrophobic slip- 
page is now widely accepted. If a liquid flows past a rough 
hydrophobic (i.e. superhydrophobic) surface, roughness 
may favor the formation of trapped gas bubbles, result- 
ing in a large slip length [9rlbJ|. For rough hydrophilic 
surfaces the situation is much less clear, and opposite 
experimental conclusions have been mad_e^ one is that 
roughness generates extremely large slip 
that it decreases the degree of slippage 



15 
If 



and one is 
13]. More 



recent experimental data suggests that the description of 
flow near rough surfaces has to be corrected, but for a 



separation, not slip [18|. The theoretical description of 



such a flow represents a difficult, nearly insurmountable, 
problem. It has been solved only approximately, and only 



for a case of the periodic roughness and far-field flow with 
a conclusion that it may be possible to approximate the 
actual surface by a smooth one with the apparent slip 
boundary condition lil- 21\. 



In this letter we address the fundamental, but still open 
questions (i) whether the effect of random roughness on 
the flow may be represented by replacing the no-slip con- 
dition on the exact boundary by an effective condition 
on the equivalent smooth surface, (ii) where this smooth 
surface is located, depending on geometric parameters of 
roughness, and (iii) does this effective condition represent 
that of slip or no-slip? We will quite generally assume 
that the flow near and far from the interface is a stable, 
laminar flow field. 




FIG. 1: (Color online) Sketch of the system: a sphere of 
radius R approaches a rough surface with a fixed area fraction 
<j> covered by roughness elements. The separation h is defined 
on top of the surface roughness at position X2 = r. 

General idea and models.— To address these is- 
sues we analyze the hydrodynamic interaction between a 
smooth sphere of radius R and a rough plane (see Fig. [J). 
Beside its significance as a geometry of SFA/ AFM dy- 
namic force experiments, this allows us to explore both 
far and near-field flows in a single "experiment" . As an 
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initial application we study roughness elements of a fixed 
height r that are distributed at random uncorrelated po- 
sitions with a given probability </>. Such a surface mimics 
a situation explored in recent experiments [l5-17|. In 
Cartesian coordinates x = (xi, X2, £3), a separation h is 
defined on top of the roughness, xi = r, which finds its 
definition in the AFM experiment [TBI, [T3] ■ 

The exact solution, valid for an arbitrary separation, 
for a sphere approaching a smooth plane is g iven by the- 



oretical solutions of Brenner and Maude 22 , 2 
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l§7 = 4sinh£ (1) 

n(n+l)[8e (2 " +1)5 +2(2n+3)(2n-l)] 
(2n-l)(2n+3)[4sinh 2 (n+±)£-(2n+l) 2 sinh 2 £] 

oo J i(Ti+l)[(2n+l)(2n-l)e 2£ -(2ra+l)(2n+3)e- 2? ]\ 
— 2^n=l (2n-l)(2n+3)[4sinh 2 (n+i)£-(2n+l) 2 sinh 2 £] J ' 

with Fst = dTtfJ-Rv, where p is the dynamic viscosity, v is 
the velocity, and cosh£ = h/R, £ < 0. The leading term 
of this expression can be evaluated as 



F2_ 9R 
F St ~ 8 h ' 



(2) 



At large separations, h 3> R, the hydrodynamic force 
on a sphere turns to the Stokes formula, but at small 
distances, h <C -R, the drag force is inversely proportional 
to the gap, F2/Fst —> 9R/(8h). A consequence of this 
lubrication effect is that the sphere would never touch the 
wall in a finite time. The flow in the vicinity of a rough 
surface should deviate from these predictions. A possible 
assumption is that the boundary condition at the plane 
X2 = r should be written as a slip condition [l5j |. To 
investigate this scenario we suggest to present a force as 
a product of Eq. [2] and a correction for slip 
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where this correction, /*, is taken to be equal as pre- 
dicted for a lubrication force between a no-slip surface 
and a surface with partial slip 24 1 . 
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Another assumption would be that the rough surface is 
hydrodynamically equivalent to a smooth one located 
somewhere between the top and bottom of rugosities (at 
X2 = r e fi = r — s). As found in (l8|). the force can be 
represented as 
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(4) 



At small h expressions [3] and |4] give different asymp- 
totic behavior of a drag force, F 3 /Fst — > 9R/(32h), and 
Fi/Fst — > 9R/(8s). While the second scenario allows a 



sphere to touch a plane, in the first model this is impos- 
sible since the drag force diverges (but differs from the 
standard lubrication asymptotics by a factor of 4). Thus, 
a drainage study allows to distinguish between these two 
models of hydrodynamic flow past rough surfaces. 

Simulation method.— We apply the lattice Boltz- 
mann (LB) method to simulate the flow field between a 
smooth sphere approaching a rough plane [H^. The 
method allows precise measurements of the force acting 
on the sphere and to explore the very large range of pa- 
rameters. Besides that, in our simulations we can con- 
sider a "clean" situation of a hydrodynamic force and 
avoid effects of surface forces which significantly compli- 
cate the analysis of SFA/AFM data. Since the method 
is well established, we only shortly describe it here. By 
using a discretized and linearized version of Boltzmann's 
equation 

ni(x + Ci,t + 1) - nj(x,i) = ^ -V,!"'/' - ^-(x,*)), (5) 
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the LB approach allows to fully resolve the hydrodynam- 
28j . Positions x are discretized on a 3D lattice with 
discrete velocities Ci pointing to neighboring sites. 
Each Ci relates to a single particle distribution function 
rii(x, i) which is advected to neighboring sites at every 
time step. Then, rij(x, t) is relaxed towards a local equi- 
librium n° q (p, j) with a rate given by the matrix elements 
Ajj. Mass p and momentum j as given by moments of 
rii(x,t) are conserved. We use the natural units of the 
system, i.e. the lattice constant <5x for the length and the 
time step St for time. Massive particles are described by 
a continuously moving boundary which is discretized on 
the lattice. Momentum from the particle to the fluid is 
transfered such that the fluid velocity at the boundary 
equals the particle's surface velocity. Since the momen- 
tum transferred from the fluid to the particle is known, 
the hydrodynamic force can be recorded. If not stated 
otherwise, the 256 3 <5x 3 system contains a sphere with ra- 
dius R = 16(5x which is moved in y direction at constant 
velocity v = 10~ 3 Sx/St. The fluid density is kept con- 
stant and the kinematic viscosity is fj,/ p = 0.1 resulting 
in a Reynolds number Re = 0.16. No-slip surfaces are 
described by mid-grid bounce back boundaries and a slip 
boundary is implemented by a repulsive mean-field force 
acting between fluid and surface @, El. We carefully 
checked the influence of system size, radius, and separa- 
tion to insure that finite size and resolution effects are 
negligible [23]. Also, by testing different resolutions we 
assured that a lateral width of roughness elements of lSx 
is sufficient. 

Results and discussion.— We test our method 
by measuring the hydrodynamic interaction between 
smooth surfaces. Fig. [2] shows the normalized hydro- 
dynamic force for two simulation sets. In the first one, 
a sphere of R = 8Sx is driven with v = 10~ 4 (5x/<5t. In 
the second run the sphere is twice as large, R = 16<5x, 
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FIG. 2: Hydrodynamic force acting on a sphere with radius 
R = 165x (triangles) and R = 8Sx (circles) driven to a smooth 
wall with velocity v — W~ 3 Sx/St and v = W~ 4 5x/8t, corre- 
spondingly. The solid and dashed curves are calculations of 
the force expected with no-slip boundary conditions at the 
wall (Eqs.[T][2]|. Squares are the results measured for a rough 
wall \cj> = 4%,r = WSx). 
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and the driving velocity is an order of magnitude larger, 
v = 10~ 3 <5x/<5t. Fig. [5] includes the exact theoretical 
curve, Eq. [TJ The fit is excellent for all separations, in- 
dicating that large shear rates do not induce any slip, a 
conclusion which does not support recent experimental 
data [|| . Note that the first-order approximation, Eq. [21 
practically coincides with the exact solution. These simu- 
lations demonstrate that finite size effects and resolution 
effects can be well controlled: for h < 2R a 256 3 <5x 3 sys- 
tem is found to be sufficient to avoid artefacts at large 
separations h (29j . Separations < lSx are excluded from 
the analysis since the finite resolution leads to larger de- 
viations. Also included in Fig. [2] is a normalized force 
measured near a rough wall = 4%, r = 10<5x), which at 
small distances is much smaller than predicted by Eq. [T] 
This is qualitatively consistent with the AFM observa- 
tions 0, [l5| , but in contrast to the SFA data [TtJ , which 
likely reflects a different way of a definition of zero sepa- 
ration in the SFA (at the bottom of asperities) . 

To examine the significance of roughness more closely, 
the force curves from Fig. [5] are reproduced in Fig. [3] 
in different coordinates. Figs. [3h. and b are intended 
to indicate that both near field and far field theoreti- 
cal asymptotics for smooth surfaces are well reproduced 
in simulations. Figs. [3}; and d show that simulation data 
for a rough surface (0 = 4%, r = lOSx) show deviations 
from the behavior predicted by Eq. [T] A possible ex- 
planation for this discrepancy is that we invoke slippage 
at the wall, as modeled by Eq. [3J This is illustrated in 
Figs. [3J: and d, where the simulation data are compared 
with another theoretical calculation in which a constant 
slip length of b = 2.55<5x, obtained from the best possi- 
ble fit of the force curve, is incorporated in the model. 
This has the effect of decreasing the force, and it pro- 
vides a reasonable fit to the data down to h/R ~ 3, 
but at smaller gap the simple model of slip fails to de- 
scribe simulation data, by predicting a larger force and 
its different asymptotic behavior. This suggests, that it 
can only be considered as a first approximation, valid at 
large distances from the wall. This conclusion is con- 



FIG. 3: Hydrodynamic force plotted in different coordinates. 
(a,b) the data sets for smooth surfaces reproduced from Fig[2] 
(c,d) show the force for two rough planes with r = 105x, 
cj> — 4% and <j> = 50%. For cj> = 4% the asymptotic behavior 
for small h cannot be fitted with a slip, but the assumption of 
an effective boundary position holds. The values for <j> = 50% 
recover the case of a flat surface at r. The data for a smooth 
slip boundary confirms the validity of Eq. [3] 



sistent with early results obtained for a far field situa- 
tion but does not support recent AFM data [15| . 
However, as shown by the simulation data, Eq. [3] is well 
applicable in the case of a slippery wall. An alternative 
explanation for the smaller force compared to the the- 
ory for smooth surfaces, Eq. [TJ can be obtained if we 
assume that the location of an equivalent effective wall, 
where no-slip boundary conditions are applied, should be 
shifted, as modeled by Eq. [4] A corresponding theoreti- 
cal calculation of the drag force is shown in Figs. [3J: and 
d. This estimate requires knowledge of the effective wall 
position r e g. The value r e g = 7.86<5x was obtained from 
the fit of the measured force curve and is enough to give 
a good match to the data at very small distances, which 



confirms the conclusions of a recent experiment [18 
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FIG. 4: Effective height r e e normalized by the maximum 
height r as a function of a density of roughness elements, 
4>, for r = 105x and r — 20<5x plotted in different scales. 

By performing similar fits for a variety of drainage runs 
with different <fi and for surfaces with different height of 
roughness elements (r — l(Wx and r = 20r5x) as well as 
its different lateral width (r5x and 2r5x) we find that the 
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same conclusion is valid for all situations, but r e g/r is it- 
self a function of <f> (being surprisingly insensitive to the 
value of r). In Fig. 3] we examine this in more detail. The 
simulation data show that r e g required to fit each run in- 
creases from to r very rapidly, so that at <j> = 20% it is 
already above 0.9r, and at <fi = 50% it is almost equal to 
r. This is illustrated by including the data obtained for a 
larger density of roughness elements {(f) = 50%, r = 105x) 
in Fig. [3b, that do not show a discernible deviation from 
the theoretical predictions for smooth surfaces. Thus, 
a small number of roughness elements has enormous in- 
fluence on film drainage, confirming earlier theoretical 
ideas [30j . 

Conclusion.— We have presented lattice Boltzmann 
simulations describing the drainage of a liquid confined 
between a smooth sphere and a randomly rough plate. 
The measured force is smaller than predicted for two 
smooth surfaces if the standard no-slip boundary condi- 
tions are used in the calculation. What our results show, 
however, is that at small separations the force is even 
weaker and shows different asymptotics than expected if 
we invoke slippage at the smooth fluid-solid interfaces. 
To explain this we use the model of a no-slip wall, lo- 
cated at an intermediate position (controlled by the den- 
sity of roughness elements) between top and bottom of 
asperities. Calculations based on this model provide an 
excellent description of the simulation data. Besides this, 
by proving a correctness of this simple model to describe 
flow past a randomly rough surface, we have suggested 
a validity of a number of simple formulas for a hydrody- 
namic drag force. Although formally they can only be 
considered as first-order approximations, their accuracy 
is confirmed by simulation. Our results open the possi- 
bility of solving quantitatively many fundamental hydro- 
dynamic problems involving randomly-rough interfaces, 
including contact angle dynamics, coagulation and more. 
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